Effects of interladder couplings in the trellis lattice 



Jose A. Riera and Sergio D. Dalosto 
Instituto de Fisica Rosario, Consejo Nacional de Investigactones Cientificas y Tecnicas, y Departamento de Fisica, 
Universidad Nacional de Rosario, Avenida Pellegrini 250, 2000-Rosario, Argentina 

(February 1, 2008) 

Strongly correlated models on coupled ladders in the presence of frustration, in particular the trellis 
lattice, are studied by numerical techniques. For the undoped case, the possibility of incommensurate 
peaks in the magnetic structure factor at low temperatures is suggested. In the doped case, our 
main conclusion for the trellis lattice is that by increasing the interladder coupling, the balance 
between the magnetic energy in the ladders and the kinetic energy in the zig-zag chains is altered 
leading eventually to the destruction of the hole pairs initially formed and localized in the ladders. 

PACS: 75.50.Ee, 71.10.Li, 75.40.Mg, 71.27.+a 



I. INTRODUCTION 

The study of strongly correlated systems in low spa- 
tial dimensions is nowadays the center of an intense effort 
both theoretically and experimentally. Among these low- 
dimensional systems the ones containing two-leg ladders 
have received considerable attention. One of the original 
motivation for the study of ladders was the search for 
a simple mechanism for pairing in a strongly correlated 
model.El This possibility was later confirmed by a num- 
ber of studies and the symmetry of the superconducting 
order parameter was found to be d^2_y2 (Refs. A 
considerable interest in these theoretical predictions was 
renewed by the discovery of superconductivity in the lad- 
der compound Sri4Cu2404i (14-24-41) after being doped 
with Ca and under a pressure of « 3GPa (Ref. |). 

However, the 14-24-41 compound, as well as many 
other compounds, like another cuprate SrCu203 and the 
vanadates CaV205 and NaV205, actually contain lay- 
ers of two-leg ladders which are coupled by frustrated 
effective interactions in the so-called trellis lattice. The 
strength of these frustrating interladder couplings may 
be weak enough to consider the ladders as essentially iso- 
lated or strong enough to change radically the physical 
behavior of a single ladder. 

The original experiments on Sro.4Cai3.6Cu2404i 
(Ref. ^ reveal that the superconducting critical tempera- 
ture reaches its maximum of 12K at a pressure of « 3GPa 
and then decreases as the pressure is further increased. 
Similar results were obtained for Sr2.5Caii.5Cu2404i 
(Ref. 1^). The application of pressure to this compound 
may change the strength of some couplings, or to addi- 
tionally increase the doping of the ladder layers as holes 
are transfered from the chain layers also present in this 
compound. The main purpose of this paper is to ana- 
lyze the first possibility, neglecting more radical changes 
in the crystallographic structure.El In a highly simplified 
model for this compound, we consider three sets of cou- 
plings: along the ladder legs, along the ruiigs, and on the 
"zig-zag" interactions between the laddersJj The effect of 
varying the rung interactions, keeping fixed the leg ones, 
on a single ladder, has been analyzed extensively.a In the 



present study we will concentrate on the effects of vary- 
ing the interladder couplings specially on the magnetic 
and pairing properties. 

To this purpose we study the t — J model which is ap- 
propriate to describe these cuprates and vanadates char- 
acterized by large on-site Coulomb repulsion and close to 
half-filling. We start our studies with undoped furstrated 
coupled ladders. For the trellis latttice, it was suggested 
that it is possible a transition from a spin liquid to a pos- 
sible spiral order with incommensurate magnetic |eorrela- 
tions as the interladder coupling (ILC) increasesB Then, 
we will consider specially the case of two-hole doping, 
i.e. the evolution of pairing as the interladder coupling 
is varied. Preformed itw)^ pairs are already present in 
the uncoupled laddersEl'cnlJ and a small ILC could lead 
to superconductivity (SC) as a proximity effect between 
the ladders. However, our main concern in this work is 
not the onset of SC but rather the effect on pairing due 
to somewhat large ILC. 

There are further motivations to study both experi- 
mentally the 14-24-41 compound and theoretically the 
t — J model on the trellis lattice. It is in effect remark- 
able how a relatively small difference like the one between 
the square lattice of the Cu-0 planes in high-Tc cuprates 
and the trellis lattice in 14-24-41 leads to such a con- 
siderable difference in the superconducting properties of 
those materials. This difference is even more_Lemarkable 
if we take into account the presence of stripesEll in the un- 
derdoped regime and (at least) above the superconduct- 
ing region of the high-Tc cuprates. These stripes can 
be thought as metallic ladders separated by insulating 
antiferromagnetic ones, specially in the "bond-centered" 
stripes-pbtained from a numerical study of the 2D t-J 
modelB. 

In this sense, the t — J model on the trellis lattice is a 
testing ground for the study of the competition between 
magnetic and kinetic energies which is at the core of the 
mechanism of micro-phase separation leading to the for- 
mation of stripes in the high-Tc cuprates.E3 Our main 
conclusion for the trellis lattice is that by increasing the 
interladder coupling, the balance between the magnetic 
energy in the ladders and the kinetic energy in the zig-zag 
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chains in between the ladders is altered leading eventu- 
ally to the destruction of the hole pairs initially formed 
and localized in the ladders. We also suggest the possi- 
bility that the hole pairs may go to the zig-zag chains in 
a process which represents a transition from a magnetic 
to a kinetic mechanism of pair binding. 

The t-J model on the trellis lattice is given by the 
Hamiltonian: 



rung 'Winter (1) 



where 



"Ha = - ^ {c\^Cj„ + h.C.) 

The couplings are (t, J), it',J') and {Unter, Jinter) along 
the legs, along the rungs and between the ladders re- 
spectively (Fig. |l|). The rest of the notation is standard. 
Periodic boundary conditions in both directions are con- 
sidered except otherwise stated. We adopt J — 0.4t, 
a value usually taken to model high-Tc cuprates, and 
in order to reduce the number of independent variables 
Ja — J{ta/t)'^- Moreover, we take t — 1. For the un- 
doped compound, njeutron scattering experiments for the 
14-24-41 compounclll3 suggest J = 2J'. However, taking 
into account the analysis mentioned in Ref . ^ most of 
our calculations have been done for the isotropic case. 

We use various numerical techniques like quantum 
Monte Carlo (QMC), with a conventional world- line al- 
gorithm, and exact diagonalization with the Lanczos al- 
gorithm (LD), complemented by the continued fraction 
formalism to compute dynamical properties. 
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FIG. 1. Left: trellis lattice. Right: modified trellis lattice. 
Ladder legs and rungs are indicated with full and dashed lines 
respectively. In the modified lattice, "perpendicular" and "di- 
agonal" interladder couplings are indicated with doted and 
dash-dotted lines respectively. 



II. FRUSTRATED COUPLED SPIN LADDERS 

The first issue we want to address is the evolution of 
the magnetic order in the absence of doping as the in- 
terladder coupling is increased. To this end we have 
computed the static magnetic structure factor S(q) on 
L X L clusters using QMC. Due to the presence of frus- 
tration, the minus sign problem prevents us to reach low 
temperatures and to study larger clusters which would 
be necessary to perform a finite size scaling. The same 
problems were already faced in aj)revious QMC study of 
the susceptibility of this system.E3 In order to reduce the 
minus sign problem and to take advantage of the simple 
checkerboard decomposition,^ we take a slightly modi- 
fied lattice in which the interladder couplings are "per- 
pendicular" (Jperp) and "diagonal" {Jdiag) as shown in 
Fig. [|. This lattice contains just half of the diagonal in- 
terladder couplings than in the trellis lattice, and we have 
checked by exact diagonalization on small clusters that 
this difference does not change qualitatively the results. 
Even for this modified lattice the minus sign problem is 
severe as shown in Fig. ||. We recall that the average of 
any observable O, in_fi system which presents this prob- 
lem is computed as:l 



{O) 



(OSign) 
{Sign) 



with respect to a modified partition function of the 2-1-1- 
dimcnsional problem, Z' = |expS'(s)| where S{s) is 
an effective action. In particular, the average sign is: 

{Sign) = -^^Sign{s)\expS{s)\, 



where Sign{s) — sign{exp S (s)) . Then {Sign) is the 
ratio of the original partition function Z = '^sGxpS{s) 
to Z' . In practice, at each measurement step, expS'(s) is 
computed as the product of the transition elements of 
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FIG. 2. Average sign of the QMC on the modified trellis 
lattice, (a) at J = J' = 1, Jperp = Jdiag = 0.2 for 4 X 4 (open 
circles), 6x6 (full circles), 8x8 (squares), 12 x 12 (diamonds) 
and 16 x 16 (stars) and (b) for the 8x8 lattice, J = J' = 1, 
Jperp = 0.4 and Jdiag = 0.4 (opcn circles), 0.8 (full circles), 
1.2 (squares), 1.6 (diamonds), 2.0 (stars). 

all tha— cubes that makes up the 2 + 1-dimensional 
lattice.EJ Nonetheless, although we cannot compute some 
quantities in the bulk limit, we can indicate qualitatively 
the behavior of the magnetic order as a function of Jmter 
(we take J as the unit of energies in this section) . 

Typical results are shown in Fig. ^ for coupled isotropic 
(J' = J) ladders on the 8x8 cluster. In Fig. ^(a) we 
show for Jperp = Jdiag = 0.2 the characteristic structure 
factor of isolated ladders with a peak at (q^, Qy) — (tt, tt) 
{x (y) is the direction along (transversal) to the ladders). 
This peak becomes more pronounced as the temperature 
is lowered. On the other hand, keeping Jperp — 0.2 and 
as Jdiag is increased, the peak starts to shift from (tt, tt) 
to (7r,7r/2) (Fig. ^(b,c,d)). A second interesting feature 
should be noticed in Fig. ^(d): the peak of S(q) is lo- 
cated at (tt, tt) at high temperature (in this case down 
to T w 0.8, in units of J) and as the temperature is 
lowered it starts to shift away from (7r,7r). At low tem- 
peratures (in this case below T « 0.4 the peak is located 
at (tt, 7r/2). Since it is clear that this behavior is caused 
by the frustration of the intcrladder couplings, it will 
certainly be present in the original trellis lattice. As in- 
dicated in Fig. |3|:, it is possible that an incommensurate 
peak across the ladder direction could be present at inter- 
mediate values of Jdiag and intermediate temperatures. 
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FIG. 3. Magnetic static structure factor S(q) obtained 
by QMC on a 8 X 8 cluster of the modified trellis lattice, at 
several temperatures (in units of J), J = J' = 1, Jp^rp = 0.2 
and (a) Jdiag = 0.2, (b) 0.4, (c) 0.6 and (d) 0.8. 

The second point we want to examine is the behavior 
of the excitations of these systems, in particular the S = 
1 excitations as can be measured by neutron scattering 
experiments. For this purpose, using conventional LD 



with the standard continued fractions formalism,t£l we 
have computed the zero temperature dynamical structure 
function {zz component) S{q,uj). In this case, we have 
to limit ourselves to somewhat smaller clusters but we 
are confident that the qualitative features we found will 
survive in the bulk limit. 

Results obtained for the 4x4 cluster are shown in 
Fig. ^. In the absence of frustration (Fig. ^(a)) the peak 
in 5(q, w) which corresponds besides to the lowest ex- 
citation, is located at (tt, tt), as expected in the bulk 
limit for an AF order. As a frustrating ILC is increased 
(Fig. ^(b,c,d)) it can be seen that considerable spectral 
weight is transferred to the peak at (7r,7r/2), which be- 
comes finally the lowest energy excitation. Similar results 
are also shown for the 4x6 cluster for q = (tt, tt) and 

The results shown in Figs. |3| and [J are unequivocally 
due to frustration and are qualitatively similar to the 
ones previously obtained ii>-& system of ferromagneti- 
cally (FM) coupled ladders.ll3 Similar results have been 
obtained by LD on the 4x4 and 4x6 clusters of the 
real trellis lattice. In the case of the trellis lattice, as in 
the FM ILC case, we expect that the behavior above dis- 
cussed will be present in the bulk limit for strong enough 
interladder couplings and low enough temperatures. The 
impossibility of assessing finite size effects prevents us to 
determine if this behavior is present for arbitrarily small 
values of Jinter Or, on the contrary, only for values larger 
than a "critical" one. 
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FIG. 4. Dynamical structure factor ^(q, uS) for several 
momentum for J), J = J' = 1, Jinter ~ 0.8 and various 
values of Jdiag as indicated in the figures. These results were 
obtained on a system of two coupled 2x4 (full lines) and 2x6 
(dashed lines) ladders. 

We have not detected awy sign of incommensurability 
along the ladder direction.EI Such an incommensurability 
is expected in principle since a trellis lattice can be also 
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considered as coupled Ji — J2 chains (Ji = Jmter, Ji = 
J, in our notation) which are known to present a peak 
in S'(q) at a momentum which continuously varies from 
TT to 7r/2 (as defined on ouunodified trellis lattice) as 
J2/J1 goes from 00 to zcro.ll3 However, notice that, as 
can be seen in Fig. ^, values of Jinter > 0.6 at low enpwgh 
temperatures cannot be reached in our simulations .c3 



III. DOPED TRELLIS LATTICE. 

We now analyze the hole pairing in the doped trellis 
lattice as the interladder couplings are increased. To gain 
some insight in this problem we start by considering an 
isolated building block of the trellis lattice. This minimal 
system is a three chain cluster consisting of a ladder and 
a zig-zag chain, (fines "1", "2" and '3' in Fig. 0). The 
Hamiltonian is the one defined in Eq. (|l|) . All the results 
in this section have been obtained by exact diagonaliza- 
tion. The justification of this study involving somewhat 
small clusters is based on an extensive body of similar 
studies of strongly correlated models which shows that 
an important part of the physics of these models is dom- 
inated by short range effects, which are appropriately 
captured in these small cluster calculations. 

In this minimal trellis lattice already appears, upon 
doping with two holes, the main feature we want to em- 
phasize. In Fig. H, the relative hole occupancy (or prob- 
ability of finding a hole) on each chain is shown in the 
3x6 cluster with two holes for t' = 0.75, 1.0 and 1.5 as a 
function of tmter- In the three cases, for small Unter the 
holes are almost completely located in the ladder legs. 
As tinter iucrcascs the probability of finding a hole in the 
outer ladder leg decreases while increases the occupancy 
of chain "3" which is connected to the ladder by the zig- 
zag interaction. There is a neat change of behavior, from 
a situation in which the occupancy of chain "3" is virtu- 
ally unoccupied to a situation in which the unoccupied 
chain is "1". moreover, this crossover is rather abrupt, 
specially for t' = 0.75. 




FIG. 5. Probability of finding a hole in the outer (circles) 
and inner (squares) ladder legs, and in the third chain (trian- 
gles) of the 3x6 cluster with two holes for (a) t' = 0.75, (b) 
t' = 1.0 and (c) t' — 1.5. The probability of finding a hole in 
chain "1" and the other hole in chain "2" , and the probability 
of finding a hole in chain "2" and the other in chain "3" are 
indicated with diamonds and stars respectively. 

For relatively small tinten the hole distribution is typ- 
ical of that of an isolated laddeiEj, i.e. they form a 
bound pair with one hole on each leg. On the other 
hand, for larger values of tinter {tinter > 0.6 for t' = 1.0 
and tinter > 0.9, for t' — 1.5, on the 3x8 cluster) the 
holes have moved from the ladder to the two chains con- 
taining the zig-zag interactions. As we indicated in the 
previous section, these two chains with the zig-zag in- 
teractions form a chain with first and second neighbor 
interactions, in this case with L^^ J couplings, i.e. a frus- 
trated ti ~ t2 — Ji — J2 chain.ES 

To understand the mechanism that produces this 
change in hole pairing, let us examine the contributions 
to the total energy from different terms of the Hamilto- 
nian (|l|) as the interladder hopping is increased. As it 
can be seen in Fig. ^, the main differences between the 
contributions from the ladder and those of the zig-zag 
chain are (i) in the ladder the magnetic energy domi- 
nates (along the legs for t' = 1.0 or along the rungs for 
t' — 1.5) while in the zig-zag the kinetic energy is the 
most important, and (ii) the main gain in energy as tinter 
is increased comes from precisely the hopping term of the 
zig-zag ILC while the magnetic energy on ladders is the 
most strongly decreased. 
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FIG. 6. Energy contributions from different terms of 
the Hamiltonian along the legs (circles), rungs (squares) and 
zig-zag (triangles). Full (empty) symbols correspond to mag- 
netic (kinetic) energies in absolute value. Results obtained on 
a 3 x 6 cluster for (a) t' = 0.75, (b) t' = 1.0 and (c) t' = 1.5. 

The gain in kinetic energy on the frustrated chain with 
respect to the ladder can be explained by qualitative ar- 
guments as is shown schematically in Fig |^. In the frus- 
trated chain we have assumed an AF order of the spins 
along the chains which is expected for J2 > Ji (^2 > ti). 
When the hole moves, as in a simple t — J chain, the 
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hole leaves behind just a single frustrated (ferromagnetic) 
bond. Something similar occurs in the case of AF order 
along the zigzag chain (Ji > J2). In the case of ladders, 
we have assumed a magnetic background formed by spin 
singlets on the rungs. As a hole moves from its initial po- 
sition, it leaves behind a string of higher energy singlets 
on the diagonals of the plaquettes. Hence there 



(a) 





(b) 



FIG. 7. Schematic picture of the movement of a hole (a) in 
a frustrated t — J chain and (b) in a f — J ladder. In each 
case, the original (final) position of the hole is shown in the 
top (bottom) panel. 

is a cost in energy which increases roughly linearly with 
the distance traveled by the hole. 

What are the possible consequences of this behavior 
found in the three-chain cluster for the trellis lattice? In 
this case, of course, any chain along the ladder direction 
belongs at the same time to a ladder and to a frustrated 
interladder chain. The question is if the holes, initially 
paired on a plaquette in isolated ladders, would tend to 
break the pairs and move independently on the frustrated 
chains as the ILC are increased. There is yet another in- 
teresting possibility that is that the holes form pairs on 
the frustrated chains. These pairs would have more ki- 
netic energy than the ones formed on ladders and this 
change of pairing would imply a change from a magnetic 
binding on ladders to a "kinetic binding" on chains. In 
any case, taking into account the results from the three- 
chain cluster, we predict a loosening of the pairing on 
ladders. It is difficult to answer these questions by cal- 
culations on finite clusters. Exact diagonalization results 
for two holes on the 4x4 cluster give support to the 
above mentioned possibilities. In Fig. || the hole-hole 
correlation functions for = 1.0 at several distances as a 
function of tinter are shown. At small tinter, the largest 
correlation corresponds to a pair of holes along the diag- 
onal of a plaquette, which is typical of isolated ladders. 
Around Unter — 0.7 there is an abrupt change to a situa- 



tion in which the largest correlations correspond to holes 
belonging to the same frustrated chain. The second hole 
is slightly more likely to be in the other chain of the same 
interladder zig-zag chain (site '4' of Fig. ||). Somewhat 
smaller is the correlation on the same chain but at the 
largest distance on this cluster (site '3'). This behavior 
is radically different to ladders AF coupled in a square 
lattice without frustration. In this case, the d-wave pair 
typical of a ladder evolves smoothly to the rotationally 
invariant d-wave pair of the square lattice. 
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FIG. 8. Hole-hole correlation functions vs. tinter on the 
4x4 cluster for t' = 1.0. The relative distances from the 
reference site are indicated in the plot. 

Similar results are obtained for the 6x4 cluster with 
two holes. In Fig.^, we show pictorially the most likely 
hole probability for t' — tinter — 1- The area of the circles 
is proportional to the probability of finding a hole if there 
is a hole in a given site. In this case the largest probability 
corresponds to holes located at the maximum distance 
along the same chain. The next probability in decreasing 
order also corresponds to a hole in the frustrated zig-zag 
chain. On this cluster we found that as tinter is increased 
from zero the holes initially at a distance ^2 starts to 
move away on the same ladder and finally they move to 
the same ILC chain. 
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FIG. 9. Most likely hole distribution if a hole is located in 
the position indicated by a solid circle. Results obtained on 
the 4x6 cluster for t' = 1.0 and tinter ~ 1.0. 

Finally, for the sake of completeness, we show in Fig. 
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the largest nearest and next nearest neighbor spin-spin 
correlations for small and large values of tinter for the 
most likely position of the holes in each case (see Fig. ||) . 
The structure of these correlations globally agrees with 
the schematic picture of Fig. M. 




FIG. 10. Largest magnetic correlations for the most likely 
hole distribution. The line thickness of each segment is pro- 
portional to the AF correlation of the connected spins. FM 
correlations are indicated with dashed lines. Results obtained 
on the 4x4 cluster for t' — 1.0, (a) tinter ~ 0.4 and (b) 

tinter — 1.0. 

IV. CONCLUSIONS 

In summary, we have performed numerical studies on 
strongly correlated electron systems, as described by the 
t — J model, on frustrated coupled ladders, in partic- 
ular on the trellis lattice. In the undoped case, QMC 
simulations, although hampered by the minus sign prob- 
lem, allowed us to reach large enough clusters to detect 
meaningful changes in the magnetic properties of this 
system. In this case, for ladders coupled with frustrating 
interactions, we have shown that the peak of the mag- 
netic structure factor shifts from (7r,7r) to (7r,7r/2) for 
low enough temperatures as the ILC is increased. More- 
over, the peak at (tt, 7r/2) becomes also the lowest energy 
excitation. This behavior is very similar to the one pre- 
viously found for FM coupled laddersll3. We have shown 
also that this behavior appears due to the onset of frus- 
tration and hence we expect that it will appear in the 
trellis lattice as well. This behavior could be detected 
experimentally on a ladder compound like SrCu203E3 or 
Sri4Cu2404i upon a suitable application of pressure. In 
fact, after the submission of this_manuscript we became 
aware of an experimental studycj on the similar ladder 
compound CaCu203. The neutron diffraction results 
reported in this manuscript indicate that the magnetic 
structure is incommensurate in the direction of the frus- 
trated interladder interaction in the plane of the ladders. 



as suggested by the present study (see also Ref. |T^) . The 
validity of the simpler FM coupled ladders model also ex- 
plains the vanishing of the spin gapcJ as due to increasing 
interladder couplings. 

Next, we have analyzed the evolution of pairing when 
the ILC are increased. In this case the physics is governed 
by short range effects and so we studied small clusters 
with exact diagonalization. Our main result is that ILC 
suppresses pairing in ladders. Results on 4 x 4 and 6x4 
clusters indicate that holes move to the chains with first 
and second neighbor interactions formed by the legs of 
neighboring ladders and the zig-zag interactions between 
them. We have identified the mechanism of this suppres- 
sion of pairing as a gain of kinetic energy of the holes by 
moving on the frustrated chains. Even for these clusters 
we have noticed important size effects which unable us to 
determine if it appears an alternative pairing of holes on 
the frustrated chains or rather the pairing is completely 
lost when ILC are large enough and holes begin to move 
independently from each other. Finally, we would like to 
stress the radically different behavior found in this case 
with respect to that found in ladders coupled in a square 
lattice, where the d^^^y^ pairing is preserved. 
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